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We point out that solar neutrino oscillations with large mixing angle as evidenced in cur¬ 
rent solar neutrino data have a strong impact on strategies for diagnosing collapse-driven 
supernova (SN) through neutrino observations. Such oscillations induce a significant defor¬ 
mation of the energy spectra of neutrinos, thereby allowing us to obtain otherwise inaccessi¬ 
ble features of SN neutrino spectra. We demonstrate that one can determine temperatures 
and luminosities of non-electron flavor neutrinos by observing from galactic SN in massive 
water Cherenkov detectors by the charged current reactions on protons. 


I. INTRODUCTION 

The historical observation of neutrinos from the supernova (SN) 1987A at the Kamiokande Q and IMB 
j| detectors had a great impact and confirmed the basic picture of stellar collapse and SN explosion ||. 
However, the number of observed events was too small to draw definite conclusions about the explosion 
mechanism or detailed properties of SN neutrinos. Hopefully, the next SN in our galaxy will be detected by 
several massive neutrino detectors, currently operating or planned. Such observations of neutrinos from a 
galactic SN should provide us a far more detailed information on properties of SN neutrinos, as well as the 
newly formed neutron star, thus bringing substantial progress to our understanding of stellar collapse. 

The discovery of atmospheric neutrino oscillations at Super-Kamiokande (SK) JlJ implies that neutrinos 
have masses and different neutrino flavors mix, and it is supported by the results from the first long- 
baseline neutrino oscillation experiment, K2K ||. Moreover, a clear evidence for solar neutrino oscillations 
into rnu/tau neutrinos has been obtained by combining the Sudbury Neutrino Observatory (SNO) charged 
current (CC) measurement Q with elastic scattering measurement at SK [frj and the most recent in situ CC 
and neutral current (NC) measurements at SNO [||. 

Neutrino oscillations add a new “complication” to the diagnostics of SN neutrinos, since it is no longer 
true that a neutrino v a leaving the neutrinosphere with definite flavor a will be detected on Earth as the 
same neutrino species. Then, the task of extracting information such as the original of neutrino spectra in 
the SN core from terrestrial observations, requires solving an “inverse problem”. 

In this letter, we point out that rather than a “complication” the neutrino oscillation provides a new 
powerful tool for probing otherwise inaccessible features of SN neutrino spectra. This holds if the mixing 
angle responsible for the solar neutrino oscillation is large || as clearly indicated by the current solar neutrino 
data. 

We show that a high statistics observation of v e through the CC reaction v e +p —> n + e + will enable us to 
extract not only the original temperature of D e but also that of SW T ) at the neutrinosphere, as well as their 
time integrated luminosities. In order to determine these parameters we employ a y 2 method to “separate” 
two different neutrino spectra at the neutrinosphere with different temperatures and luminosities. Hereafter 
we use a collective notation v x (p x ) for v^t) {vp,(r)) because they cannot be distinguished by their physical 
properties inside the SN. 
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We focus on v e observation for the following reasons: (i) this is expected to be the channel with highest 
statistics in a large water Cerenkov detectors such as SK, Hyper-Kamiokande (HK) or UNO detectors, which 
are under consideration s Neutrinos from a galactic SN at 10 kpc would produce 7000-10000 events at 
SK and 2-3 xlO 5 events at HK [|FiJ. (ii) the large solar neutrino mixing angle necessarily implies that the 
v e spectrum observed at the Earth is a strong mixture of two originally different SN neutrino spectra of 
v e and v x . This fact has been used to derive various constraints on neutrino mixing parameters p~j| . 
Complementary information may be obtained also through direct detection of v x through the NC reaction 
either at SNO or KamLAND [ jU5| . 

In addition to generic features of SN neutrino spectra obtained in SN simulations |||l6| [jl| we will take 
into account some new features indicated by recent studies. Most importantly, they include a new parameter 
which characterize the departure from the equipartition of integrated luminosities to electron and other flavor 
neutrinos, the quantity of greatest uncertainty in SN simulations. (See below.) 


II. SUPERNOVA NEUTRINOS; BASIC PROPERTIES 


We now briefly summarize the features of neutrino spectra relevant for our work. In a SN explosion driven 
by gravitational collapse, about 99 % of the total binding energy of the neutron star, Eb ~ 3 x 10 53 erg, 
is released in the form of neutrinos during the first ~ 10 seconds after the onset of core collapse. It is 
well known that the time-dependent energy spectrum of each neutrino species can be approximated by a 
“pinched” Fermi-Dirac (FD) distribution (u|. In this work, we assume that the time-integrated spectra can 
also be well approximated by the pinched FD distributions with an effective degeneracy parameter rj 
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where E denotes the neutrino energy, and T the effective temperature. We have checked the validity of this 
approximation by using results of 20 M 0 simulation by Mayle and Wilson we take in our analysis 

rj = 2.6 and rj = 0 for D e and v x , respectively, which reproduce well the time-integrated spectra. 

Because of the hierarchy in strength of interactions with the surrounding matter neutrino temperatures 
obey (T„ e ) < (T 5e ) < (T Ux ) ~ (Tp x ) as confirmed by various SN simulations. Typical values of the average 
energies of the time-integrated neutrino spectra obtained are (E Ve ) ~ 12 MeV, (Ep e ) ~ 15 MeV and {E V;r ) ~ 
(Ep x ) ~ 24 MeV. We introduce, for later use, a new parameter te which is defined as the ratio of the average 
energies of the time-integrated neutrino spectra, 


te = 


(E9.y 


( 2 ) 


SN simulations indicate te ~ 1.25 — 2.0. 

Recent studies and SN simulations indicate various new features which were not taken into account in 
previous studies. There are three effects at least. The first is the possibility of a gross violation of equality 
in integrated luminosities of v e jv e and v x /v x by up to ~ 50 % |2l]]. The second is a violation of equality of 
physical properties of v x to v x due to the effects of weak magnetism J2^] . The third is possible difference in 
integrated luminosity between v e and v e pT| . p0p~i| | . 

The actual situation regarding violation of equipartition can be less dramatic since the ^ e i7 e annihilation 
process enhances the v x /v x luminosity, as recently shown by Buras et al. j2f|. However, we prefer to introduce 
a free fit parameter £ defined as a ratio of integrated luminosities of v x to D e \ 
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where El°J = f L Pa dt and £ = 1 in the equipartition limit. It quantifies the departure from equipartition 
of integrated luminosities to v e and v x . On the other hand, we ignore the latter two effects in this paper. 
We feel that taking into account the uncertainly in £ gives us a reasonable framework at least as a first 
approximation; the other effects are not very sizable, < 10 %, and may cancel with each other, e.g., effects 
of j2)| and Q on temperatures of v x to v x . Therefore, we treat E^, T Se , = E^. and T Vx = T Px as 
free fit parameters. 
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III. ANALYSIS METHOD 


Taking the best-fit values of mixing parameters from the latest analysis of the solar neutrino data [^5| , we 
compute the conversion probability for v e «-» v x oscillations in the SN envelope, following the prescription 
in Ref. |2(| , assuming an ap pro ximate power-law density profile, p(r) ~ r -3 . For simplicity, we neglect 
possible Earth matter effects [|27j which will depend on the location of the detectors. In the context of three- 
neutrino flavor mixing, our treatment applies to the normal mass ordering (Am 2 tm > 0) and the inverted 
mass ordering (A m 2 atm < 0) with non-adiabatic high-density resonance |l||,|28|], where A?n 2 tm denotes the 
atmospheric neutrino mass squared difference. Further implications of our proposal for the case of generic 
three flavor mixing will be given in Ref. J 20 ) . 

We take, for definiteness, the SK and the HK detectors and assume their fiducial volumes as 32 kton and 1 
Mton, respectively. We assume v e 's will come from a galactic SN located at 10 kpc from the Earth. We only 
consider v e CC reaction on proton and neglect contribution from v a e elastic scattering and CC reactions 
on oxygen in our analysis. It should give a good approximation because these processes have very small 
cross sections. The effect of weak magnetism is taken into account |d(]]. Without knowing the detection 
efficiency and energy resolution expected at HK, we assume that they are the same as those in SK and set 
the threshold energy to 5 MeV psj . 

Lacking real galactic SN neutrino data at hand, we generate an artificial data set by adopting the model 
spectra as described before. We first define arbitrarily a set of initial values for the four relevant astrophysical 
SN parameters a 0 = {E°, (Ep e )°, £°} as “true values” given by nature. The data thus generated, 

N° hs = Ni(a° ), are taken as the “observed” values. To quantify how well a single FD distribution can be 
discriminated from a superposition of two FD distributions we employ a y 2 minimization fit in the space of 
four astrophysical parameters a = {Eb, (E„ e ), te , £}■ The y 2 function is defined as 


N bin 

y 2 = 2 x Y, i N i( a ) - N i( a °) +N i (a°)ln[N i (a 0 )/N i (a)]} , (4) 

i=1 

where Wnn = 20. For simplicity, we present below the results obtained when event-by-event fluctuations of 
the artificial data set are neglected. However, we have explicitly verified that our results agree well with the 
ones obtained by generating 10,000 sets of Gaussian fluctuated data. 

Current solar neutrino data strongly favor the large mixing angle (LMA) MSW solution. Therefore, we 
focus on the LMA solution. It also gives a conservative estimate of the oscillation effect compared to the 
LOW and the VAC solutions because the effect is larger for larger mixing angles. 


IV. RESULTS AND DISCUSSIONS 

We show in Fig. the 3 er C.L. allowed parameter region obtained in the ( Ep e ) — te and (E„ e ) — Eb planes 
for SK and HK detectors by assuming the LMA solution, taking the best-fit value for mixing angle, tan 2 6 = 
0.42 ^5|. We also present results for the no-oscillation case for comparison. We set the initial astrophysical 
parameters as ( Ep e )° = 15 MeV, = 1.4, E$ = 3 x 10 53 erg, and £° = 0.5. 

Fig. 0 demonstrates that we can extract the both v e and v x temperatures in the presence of large mixing 
oscillations. The accuracies we can achieve with the LMA best-fit parameters are A te/te ~ 9 (1.5) %, 
A (Ep e )/(Ep e ) ~ 4 (1) In contrast, in the absence of oscillation there is no sensitivity in Eb due to the 
inability of determining v x and v x luminosities. The improvement in the accuracy of the determination of 
these quantities for the oscillation case is remarkable, especially with HK. 

Let us now examine how accurately we can determine the equipartition-violation parameter £ and its 
correlation with te and Eb . This is shown in Fig. ||. We note that the accuracy of te determination is 
remarkably good in spite of the rather poor accuracy in the knowledge of £. It should be noticed that Fig. 
demonstrates that HK can do a much better job for £, A£/£ ~ 10 %. We emphasize that without having 
sensitivity to £ we can not determine Eb in a good accuracy, as they are strongly correlated (See the right 
panel of Fig. g|) We finally note that HK’s enormous sensitivity ~ 10 % may mean one could potentially 
examine the problem by treating accreting and thermal phase separately under the present approximations. 
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FIG. 1. Extracting the astrophysical parameters. The figure shows 3 a contours assuming (Ep e }° = 15 MeV, t° e 
= 1.4, E ° = 3 x 10 53 erg, and = 0.5 for the SK and HI< detectors assuming the LMA solution to the solar neutrino 
problem. Best fits are indicated by the stars. 


The negative correlation between £ and te can be understood due to the fact that the effect of lowering 
which implies a decrease of the relative contribution of v x , can be compensated by increasing te to keep 
the higher energy tail of the observed spectrum similar. On the other hand, the strong positive correlation 
between £ and Ef, just reflects the relationship Eb = 2(1 + 2£)£ , O t , which is valid under our approximation 
El°* = Ey A . The correlation is robust and exists with and without oscillations as indicated in Fig. 2. 



\ \ 


FIG. 2. Determination of the non-equipartition parameter £ and its correlations with te and Eb . Same assump¬ 
tions as in Fig. jl|. 
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We have also examined the case of a smaller value of the initial v x energy, = 1.25 to verify the robustness 
of our results. We found that although the sensitivity to E\, determination worsens by a factor of ~ 2, we 
still have a reasonable sensitivity to te, Ate/te ~ 14 % (2.5 % for HK), excluding the case of te = 1 at 3 
a. See Ref. @ for details. 

We have verified that our results do not change very much even if we relax our assumptions by taking 
into account deviations from power-law density profile of progenitor star, or from the precise best-fit values 
of solar neutrino mixing parameters we have adopted. The inclusion of possible Earth matter effects for a 
given experiment will imply some regeneration effect. Although this will somewhat weaken our results, the 
effect is rather small in practice. Additional details related to the robustness of our method can be found in 
Ref. [|§. 

In summary, we have suggested a simple but powerful way of extracting separately the temperatures of 
both v e and z/ M ( r ) as well as their integrated luminosities by analyzing v e events that would be recorded 
by massive water Cherenkov detectors in the event of a galactic supernova explosion. In particular, an 
extraordinary power of megaton-class detectors (Hyper-Kamiokande or UNO) are noticed with regard to 
the determination of the integrated luminosities of v e and ^ M ( r ) as independent fit parameters. We stress 
that the large mixing between v e and i which is clearly indicated by the current solar neutrino data, 
is essential in determining temperature as well as luminosity of iv r ), and this must have more profound 
implications which await further investigation. 


Note added: 

When the first version of this paper was ready for submission to the electronic archive we became aware 
of the paper by Barger et al. Q who pursued the similar strategy as ours. However, they do not treat the 
violation of equipartition of energies as a fit parameter. 
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